Calculating the local-type /nl f° r slow-roll inflation with a non-vacuum initial state 
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Single-field slow-roll inflation with a non-vacuum initial state has an enhanced bispectrum in the 
local limit. We numerically calculate the local- type /nl signal in the CMB that would be measured 
for such models (including the full transfer function and 2D projection). The nature of the result 
depends on several parameters, including the occupation number Nk, the phase angle 9k between 
the Bogoliubov parameters, and the slow-roll parameter e. In the most conservative case, where 
one takes 9k ~ rjok (justified by physical reasons discussed within) and e < 0.01, we find that 
< /nl < 1.52(e/0.01), which is likely too small to be detected in the CMB. However, if one is 
willing to allow a constant value for the phase angle 9k and Nk = 0(1), /nl can be much larger 
and/or negative (depending on the choice of 6k), e.g. /nl ~ 28(e/0.01) or — 6.4(e/0.01); depending 
on e, these scenarios could be detected by Planck or a future satellite. While we show that these 
results are not actually a violation of the single-field consistency relation, they do produce a value 
for /nl that is considerably larger than that usually predicted from single-field inflation. 



I. INTRODUCTION 

Most cosmologists share the belief that the universe 
underwent an early inflationary period, where it ex- 
perienced an accelerated expansion 1] that produced 
the large-scale perturbations that seeded the cosmic mi- 
crowave background (CMB) anisotropy and large-scale 
structure However, understanding the mechanism 
and exact history of this early period is still the subject 
of active research and the impetus behind the develop- 
ment of a variety of inflationary models Q . One of the 
most promising avenues for discrimination between dif- 
ferent models is the measurement of the non-Gaussianity 
of the primordial curvature perturbation £, that is, the 
degree to which the primordial fluctuations contain sta- 
tistical information beyond the power spectrum [4] . As it 
stands, all measurements are consistent with a Gaussian 
curvature perturbation however, this analysis may 
change with the release of data from the Planck satellite, 
as well as large-scale structure observations. 

The bispectrum is the lowest-order non-Gaussianity 
and, at least in most models, the most likely to be de- 
tected. It is defined in terms of the Fourier transform of 
the three-point function 

(Ck x Ck 2 Ck 3 ) = (2ir) 3 5W(j:^)B c (k u k 2 ,h), (1) 

where the ^-function is required by translation invariancc 
of the underlying theory. The bispectrum B^(ki, k 2 , k 3 ) 
is a real- valued function of three variables; however, the 
measured data has only limited signal strength so the 
bispectrum is often reported in terms of /nl parame- 
ters, each of which is obtained by fitting the measure- 
ments to a template bispectrum. This paper will deal 
with the local bispectrum J^l (henceforth simply called 



/nl), the most commonly reported non-Gaussianity pa- 
rameter. The local template is the result of assuming 
that the non-Gaussianity is produced from a background 
Gaussian field £ 5 by a purely- local process [6[ , so that 

C(x)=C 9 (x) + (3/5)/nlC 5 2 (x) + ... ; (2) 
a straightforward calculation then yields the bispectrum 

Bf c (k 1 ,k 2 ,k 3 ) = ^/NL[P C ( fc i)^( fc 2) + Pdk2)P C (k 3 )+ 



Pf(fc 3 )P C (fci)] 



(3) 
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Noting that the measured power spectrum P((k) oc fc~ 3 
and that the wavenumbers must obey k! + k 2 + k3 = 
0, we see that this bispectrum peaks when one of the 
wavenumbers is much smaller than the other two, e.g. 
k 3 <C k\ fo; this limit is known as the squeezed-limit 
or local-limit [7[. In this regime, the bispectrum becomes 

4 oc (fc 1 ,fc 2 ,fc 3 ) fc3 « fcl ^ 2 « y/MLiWiW- (4) 

/nl gained much importance when Creminelli and Zal- 
darriaga [1] demonstrated a consistency relation that 
shows that for all single-field inflation models (regardless 
of kinetic term, vacuum state, etc.) in the squeezed-limit 

B c (k 3 « fr) = (1 - n.)P(ki)P(k 3 ). 

Since the local template is most sensitive to squeezed- 
limit wavenumber configurations, the consistency rela- 
tion implies that 

/ NL «(5/12)(l-n s )«0.01, 

where n s is the spectral tilt; thus, a larger detection of 
/nl would strongly disfavor single-field inflation. The 
consistency relation has been explored in many papers, 
including [9Hl3l|. However, one must be cautious when 
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using the consistency relation to interpret observations. 
The relation is inviolable only in the exact squeezed limit 
while measurements are made over a finite range of k 
values and involve a best-fit of the measured bispectrum 
to the local-template bispectrum. Recently, Agullo and 
Parker [3] showed that there can be an enhancement 
of the local-limit bispectrum for slow-roll inflation given 
a non-vacuum initial state. In § [Til we reproduce their 
result, though in a formalism where the non-standard 
initial state is recorded in the Bogoliubov parameters 
(rather than by a density matrix). We then proceed, 
in § |TTT] to calculate the value of /nl that such a model 
would produce in the CMB. In § IIVI we report the re- 
sults, which we interpret further in § [V] We find that 
/nl S 3 (i-e. undetectable) in the most conservative case 
but that we can have /nl 3> 1 if we allow some freedom 
in the parameters; this result highlights the importance 
of properly applying the consistency relation. 



II. SLOW-ROLL INFLATION WITH A 
NON-BUNCH-DAVIES INITIAL STATE 

First, we must calculate the bispectrum for slow-roll in- 
flation with a non- vacuum initial state. The fluctuations 
in the scalar curvature perturbation £ are produced by 
quantum fluctuations in a scalar field 4> which has the 
action 



b ~ 2 



d 4 x^/=g[R - (V0) 2 - 2V {<!>)) 



We suppose that the slow-roll parameters are small, i.e. 
1 -,. 



2 /H 2 < 1 



VSR 



-0/ 110+ (1/2) <j> 2 /H 2 < 1 



After switching the degree-of-freedom to C, we can calcu- 
late the action to third order (note that M~ 2 = 8nG = 1) 
as 



a(d() 2 



(5) 



where as usual a is the scale factor and H = a/ a; note 
that we have neglected a field redefinition which is a small 
effect and irrelevant for this discussion. Then, to lowest 
order in slow-roll, the equation of motion for £ is 



(6) 



Using the in-in formalism, the tree-level three-point 
function for ( is [IH 

<c 3 (n) - -i f ^'(o|[c 3 (n,iz/(i')]|o) , (?) 



where the brackets on the LHS denote the ensemble aver- 
age, Hi is the interaction Hamiltonian (and equals — L3), 
and to is the time where we set initial conditions on the 
state of the universe. To evaluate the RHS, we quantize 
£ using Quantum Field Theory in curved space-time [l6| : 



Ck(i) = M/c(i)ak + u%(t)a[ 



k ■ 



(8) 



where ak, are annihilation, creation operators, respec- 
tively, for £ modes and u k (t) are the mode functions for £, 
i.e. complex solutions to the equations of motion (jSJ) sat- 
isfying some normalization condition. We usually choose 
u to correspond to the so-called Bunch-Davies (BD) vac- 
uum, where it contains only positive frequency modes 
and where the vacuum state is unchanged under trans- 
formations that leave the metric de-Sitter; for slow-roll 
inflation, the BD vacuum state is 



H 2 1 



= (1 + ikr,)e- lkv . 



(9) 



More recently, people have realized that the assump- 
tion of a BD vacuum may be too restrictive and have 
considered initial states containing particles [ll|, [2 [13- 
[2lj . These states can be thought of as parameterizing 
either early-time or high-energy physics and can be rep- 
resented by performing a Bogoliubov transformation on 
the BD u: 



u k (t) = a k u k (t) + fa u* k (t) 
H 2 1 



a k - 



(1 + ikr])e- lkri + 



H 2 1 



= (1 -ikj])e lk ^ 



(10) 



for slow-roll inflation. If u k is a vacuum mode function^ 
then u k corresponds to a mode containing N k = \j3 k \ 
particles in mode fe; also, the Bogoliubov parameters 
obey |afc| 2 — |/3fc| 2 = 1. Only the relative phase 9 k be- 
tween a k and j3 k is physically relevant so we can write all 
quantities involving the Bogoliubov parameters in terms 
of N k andflfc (e.g. a k (3* k +a* k /3 k = 2^/N k {N k + T)cos0 fc ). 
8 k is an additional unconstrained parameter and for very 
particular choices it can dramatically enhance /nl; this 
will be further discussed later. 

To use the non-BD u, we simply replace u by u in the 
quantization of Q {5} . For the power spectrum, this gives 
us 

P ( k ) = \Uk\l-+ = ~U ^3 \ ak + ^'l 2 

= (l + 2N k + 2^/N k {N k + 1) cos0 fe 

(11) 

While we currently can only conjecture about the exact 
form of Nk, we can place on it certain constraints: 
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(a) The energy content in the fluctuations must be finite. 
This requires that, in the UV, N k = 0{l/k 4+s ) [3 
(heuristically at least, this follows from demanding 
that the energy density J d 3 kkNk converge). Ac- 
tually, this constraint has only a limited impact on 
our results; we can only observe a restricted range of 
wavenumbers so that we are free to posit any form 
for Nk for high enough wavenumber. 

(b) A more restrictive condition comes from demanding 
that the energy content in the fluctuations be less 
than that in the unperturbed background field; oth- 
erwise, there would be back reaction on the infla- 
ton's dynamics and inflation would not proceed as 
expected. Roughly, we can write this condition as 

J cPkkNk < M^H 2 . 

(c) The amplitude of the perturbations must be compat- 
ible with the observed primordial power spectrum. 
This constraint is fairly weak since, as we can see 
from (JTTJ), P(k) cx {H 4 /<p 2 )(l + 2N k ) and we can- 
not independently constrain either the first or second 
term. 

(d) The running of Nk must conform to the observed 
spectral tilt 1 — n s s=s 0.032 [5|. This implies that 

l-n s ee -dlog(fc 3 P(fc))/dlogfc 

= -2?7sr + 6e - dlog(l + 2N k )/dlog k 
= 0.032. 

Assuming e and tjsr are small, then Nk must change 
slowly 



As in [l9| . we will consider Nk with the form 



N k « A M e- fe2/fc - , 



with some cutoff scale fc cut . Condition (jbj gives Nk.o < 
M^H 2 /k 4 ut . If we take fc cut ~ ^M p iH, i.e. the scale of 
inflation, Nkfi can be of order one; for other scenarios, 
iVfc,o might be higher or lower. A more careful calculation 
gives a comparable result [l8[ . If we take the cutoff scale 
fc cu t to be sufficiently larger than the primordial modes 
observable today then, for observational purposes, Nk ~ 
Nk.o — const; this is one case we consider. We will also 
calculate the effect of having the cutoff fc cu t lie within the 
observable modes. For some cutoffs, this situation may 
not satisfy condition fd} above but it provides an idea 
about the effect of a varying Nk- 

The form of 9k is constrained only through the power 
spectrum allowing us considerable freedom in its 

choice. If we think of the non-vacuum initial state as pa- 
rameterizing unknown physics, then we expect that it is 
determined by matching the mode function (|TU)) to the 
initial conditions. It follows that the main (or at least 
a large) contributor to 9k will be the exponential factors 
e -ikr la ^ gjfeijo evaluated at some initial time rj ; the alter- 
native is that these factors are coincidentally cancelled 
off in the initial conditions which (since the initial con- 
ditions probably represent very different physics) seems 
unlikely. Thus, it makes sense to consider 9k ~ kr/o, 
though we will also explore the effect of a constant 9k for 
very large fc cut . 



Returning to the calculation of the bispectrum, we plug ([5]) and ([5} into (J7J) and find that 

4> 4 ( l \ l * * * 

B c (fa,fa,fa) = ZijjeYl [tfj UkiityukiityUkaif)) J dfj—u'l^'lyl^c.c. (12) 
We can now replace u by u to account for our initial state, so the integral becomes 

Jr l0 V d V 8 J Vo 

(13) 



For the BD case (a = 1, (3—0), the integral 
gives [l/i(ki + fc 2 + k^e^+^+^l^, which in the 

squeezed limit becomes (l/2iki)e 2lk:iV \ r> . If we consider 
a non-BD initial state, j3k ^ 0, and we get terms of a dif- 
ferent character; for example, there is a 0^^% a % 3 term 
which gives [l/i{kx - fa + fc 3 )]e i (' Cl "' C2+ ' £3 )' ) |^, or in the 
squeezed limit fc3)e lfe3?7 |' 7 . Thus, terms like this lat- 



ter term are enhanced by a factor of k\ /fa 3> 1 in the 
squeezed limit. 

In principle, it is not clear what happens at the lower 
limit of integration, when 77 = rjQ. We expect the initial 
condition for the modes to be set while the modes are 
within the horizon, i.e. krjo 3> 1. Thus, depending on 
the exact value of fc and 770, the exponentials will oscil- 
late; when performing our calculations, we will take the 
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average value, i.e. that the exponentials are zero. For the left endpoint, as usual, rj — s- 0. 
I 



In this way, we arrive at an expression for the bispectrum 



B c (k 1 ,k 2 ,k 3 ) 



i H 6 



1 



\ " 



"4 j>2 k x k 2 k 3 ^ \kf 



x / dr, «e*" + Z^e"^) «e*" + «e ife3 " + Pl 3 e'^) + c.c. , (14) 

J no 



which is written out explicitly in the Appendix. The non- 
BD part of the shape is primarily squeezed with some en- 
hancement in the folded regime (i.e k\ « k 2 ~ (l/2)fc 3 ), 
coming from the terms mentioned above. As a check on 
our calculation, note that we recover eq. (4.4) of Mal- 
dacena's paper Q if we specialize to the BD case, where 
a k 1. U. 0. 

Thus, the non-BD bispectrum is enhanced in the 
squeezed limit, resulting in a larger /nl- 



Comparison with previous work 

This result is essentially equivalent to eq. (5.9) of Ag- 
ullo and Parker [lij]. However, Agullo and Parker used 
a density-matrix formalism where the results are essen- 
tially expressed in terms the expectation value of prod- 
ucts of the occupation number operator iVkfe. By con- 
trast, our results are expressed in terms of the Bogoliubov 
parameters ak, Pk- While the formalisms are equivalent 
in principle, the relationship between them is not trivial, 
since any non-Bunch-Davies state in terms of Bogoliubov 
parameters is an infinite sum of states with definite occu- 
pation numbers [22]. We have not encountered an explicit 
formula connecting the two descriptions. 

We feel our result is worthwhile because it connects 
to the formalism of much of the previous literature on 
the subject. Also, it is often easier to use when try- 
ing to hypothesize possible initial states since these are 
often naturally expressed in Bogoliubov parameters, as 
we saw in our arguments about the expected behavior 
of the phase angle 9k- We do however admit that Ag- 
ullo's approach has the potential advantage of being able 
to calculate the bispectrum given any arbitrary initial 
state, since the space of states expressible in terms of 
Bogoliubov parameters does not fill the Fock space. 

A similar situation was treated in [19] but they were 
expecting a signal only in the folded limit and did not 
look for one in the squeezed limit. [Note that there is 
an error in their eq. (3.16) which further obfuscates the 



squeezed limit enhancement]. [20j used a similar setup 
as the present case, examining slow-roll inflation with a 
modified initial state. Like in the current paper, they 
performed a 2D projection, though in the flat-sky ap- 
proximation, which works accurately only for I > 150 
|23j and so is not ideal for calculating the local template 
(which has a strong signal when one of the I's is small) ; 
however, they used the same incorrect bispectrum as [19j . 
and so their result missed the local limit enhancement. A 
later paper including two of the same authors [21[ found 
that models with a small speed of sound and a modi- 
fied initial state can have an enhanced local signal; they 
also noted that the size of the enhancement depends on 
the phase angle and they studied the effect of different 
(constant) phase angles. However, they did not perform 
the 2D projection so they could not calculate what /nl 
would be. Furthermore, the model they used was more 
complicated, and we feel that are advantages to studying 
this unexpected behavior in the simple case of standard 
slow-roll inflation. 



III. THE OBSERVED VALUE OF / NL 

Determining the /nl signal in the CMB requires some 
additional work. In particular, we must use a trans- 
fer function to calculate the present-day temperature 
anisotropy from the primordial curvature perturbation 
and then project the 3D CMB onto the 2D projection 
that we actually observe; this gives us the angular bis- 
pectrum. Then, we use the optimal estimator to extract 
the value of /nl; essentially, this last step involves fit- 
ting the measured angular bispectrum to the template 
angular bispectrum derived from ([3]). 

The calculations required arc described in [24]. The 
angular power spectrum and bispectrum are calculated 
from their primordial counterparts by 

Q= {^L) j dkk2p dk)9 2 Ti(k), (15) 

and 
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J hhh 



-^j J r 2 dr J dk 1 dk 2 dk 3 (k 1 k 2 k 3 ) 2 B c (k l7 k 2 ,k 3 )g T i 1 (k 1 )g T i. 



(ki)9Th (h)jh (far)ji 2 (k 2 r)j h (k 3 r) , (16) 



where gri and ji are the radiation transfer function and 
spherical Bessel functions, respectively. The angular bis- 
pectrum 6z 1 i 2 / 3 for the non-vacuum initial state is calcu- 
lated using (|14p while bfif 2l for the local template bis- 
pectrum is calculated using © . To calculate C; , one can 
use (fTTj) . Finally, the optimal estimator for /nl is given 
by 



/NL 



— V 

h>h>h>2 



I 2 



b \°hl 3 b hhh 



(17) 



where 



N 



loc 



E 



i 2 



h>l 2 >l 3 >2 

for convenience, we have defined 



&hizi 3 CuCuCu 



if all l's the same 
if 2 l's the same 
if all l's different 



L hhh 



(2/i + l)(2/ 2 + l)(2/ 3 + l) 



-in 







where the matrix above is a Wigner 3-j symbol. 

Given the complicated form of gri and of the equa- 
tions above, analytically solving for /nl is generally an 
intractable problem and must be done numerically. Even 
so, the system is somewhat involved and, in general, 
presents a difficult computational challenge. This situ- 
ation is improved if the bispectrum is separable, i.e. if it 
can be written in the form 

B c (k 1 ,k 2 ,k 3 ) = fi{ki)gi(k 2 )hj(k 3 ). 

i 

While the local form bispectrum ([3]) is of this form, the 
non-vacuum slow-roll bispectrum (|14[) is not. However, 
since calculating /nl requires fitting a bispectrum to 
the local form, the result is primarily dependent on the 
regime where the local form peeks, i.e. the squeezed limit 
^3 <C fci « k 2 . More formally, one may examine (|16p. 
(TT7)) . - noting that gTi(k) peaks when k « l/r-^s, where 
tls is the comoving distance to the surface of last scatter- 
ing - and arrive at the same conclusion. In other words, 
we can make our bispectrum (|14p separable by assuming 
that it is evaluated in the local limit, where k 3 <C k\. 
Then, our bispectrum becomes: 



B c (ki, k 2 , k 3 ) 



1 



ft 4k x k 2 k\ 



(18) 



where 



- (a «-> /3) + c.c. = 
= 27V! + 2N 2 + 4N X N 2 



+ 4N 1 ^N 2 (N 2 + 


1) cos( 


h+ 


+ 4N 2 y / N 1 (N 1 + 


1) COS ( 


h+ 


+ 2^N 1 (N 1 + 1) 


cos 6*i - 


h 2V 'N 2 (N 2 + 1) cos 


+ A y /N 1 N 2 {N 1 + 


l)(N 2 


+ 1) COS 01 COS 02 + 


+ 4y/N 1 N 2 (Ni + 


1){N 2 


+ 1) sin 0i sin 02 + 


+ 4^N 1 N 3 (N 1 + 


1)(N 3 


+ 1) sin 0i sin 03+ 



+ - 



where Ni = , Bi = 0fc 4 . (To ascertain the reasonable- 
ness of this approximation, we performed the calculation 
to the next order in k 3 and found results that were 100 
times smaller). Note that the standard 1 — n s terms are 
absent because they are lower order in k 3 . 

To find /nl, we plugged (fl8|) into (fT6|) and (fT7|) and nu- 
merically performed the integrals in (I16[) . To determine 
the radiation transfer function gTi, we used the gtFAST 
software [H], based on CMBFAST 4.0 [26]; for input, 
we assumed h = 70.3, T CM b = 2.725 K, n b = 0.0451, 
n c = 0.226, n A = 0.729, Vl v =Q,N V = 3.04, r = 0.088, 
and Zrckm = 10.4 (from [5j). We summed Vs up to 2000, 
to reflect the signal available in the CMB temperature 
power spectrum. 



IV. RESULTS 

As described in § HH we calculated /nl for various 
values of the cutoff fc C ut: 

a) k cut is very large, so that, Nt ~ Nk.o ~ const: 

We supposed either that k ~ const or 6 k ~ ki] . For 
the latter case, the terms in (Til?)) with one trigonometric 
function will clearly average out to zero in (|16|) . The 
terms with two trigonometric functions initially appear 
as if they might contribute, since k\ ~ k 2 , giving e.g. 
cos 0i cos 02 ~ cos 2 0i. However, the integrand in (TTBl) is 
sensitive only to k values that satisfy a triangle inequality 
(while this is not obvious from (T17J1) . it is clear from the 
derivation of that equation) so that, for example, k 2 £ 
[k 1 — k 3 ,ki + k 3 ] . Even though k 3 <C k% ~ k 2 , r/ is 
sufficiently large (k 3 rjo ^> because we suppose that the 
initial conditions are set when k 3 is within the horizon) 
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that these terms will oscillate greatly during the integral 
over fc 2 and average out. Thus, we will be left with only 
the first line of the second equality in fp~9|) : Me h ^kr] = 

For either form of 9j~, the numerical integral is the 
same. Performing the integral, we found 



f 7fi ^ N 1 .0 6 N 



(20) 



where we have defined P(k) = H A /(2fc 3 </r) P, so that 
P(k) equals the term in parenthesis in (fTT|) . We have 
expressed the answer in terms of the slow-roll parameter 
e. If the single-field relation r = 16e still holds for this 
situation (where r is the tensor-scalar ratio), the bound 
r < 0.2 @ yields e < 0.013; however, since we have 
a non-vacuum initial state, this relation may no longer 
hold, and e could be larger than otherwise expected (but 
the slow-roll approximation breaks down for e > 0.1). 
For the case that 9 k ~ krj , (|2T)|) becomes 



*■« 20 



FIG. 1. For very large fc cut , where N k ~ Nk,o = const, 
and fine-tuned values of 8 k = const, /nl can be large. The 
maximum values (upper solid line) are achieved for cos# = 
(1 - 127V - 12N 2 )/6^/N(N + 1)(1 + 2JV), and the minimum 
values (lower dashed line) are achieved for cos# = — 1. 



/nl - 1.52 



0.01 



1 - 



1 



l + 4iV M +4iV| 



(21) 



We thus have /nl ^ 1.6(e/0.01); for the standard bound 
on e this is indistinguishable in the CMB - since the best 
the CMB can do is A/nl ~ 3 [27]. For larger values 
of e, this might actually be observable by future CMB 
satellites. 

If we instead consider that 9k « const, we can get 
significantly larger results and even negative /nl ■ Indeed, 
finding the extrema of (|20|) with respect to a constant 9k 
for Nk o > 0.017, one finds a maximum of 



/NL 



0.01 V4 



1.52 -4- (- + 9N k n + 9N? n | (22 i 













/cut 


fccut 


N k ,o 


/NL (ok 


f 


very larj; 


;e very large 


N k ,o 


1.52[1 - (l + 47V fc . 


+ 4/VI.0)- 1 ] 






0.1 


0.14 




1500 


0.104 Mpc" 1 


1 


0.44 








10 


0.46 








0.1 


5.4 x 10" 


-3 


500 


0.0347 Mpc" 1 


1 


1.9 x 10" 


-2 






10 


2.32 x 10 


-2 



TABLE I. The observed /nl for N k = /Vfc, e" fc /k ^. The top 
line shows the result for fc cut very large, so that Nk ~ const. 
For cutoff fc cu t = icut/^LS, we get a suppression of Nk for the 
wavenumbers that contribute primarily to multipoles larger 
than /cut in the CMB temperature anisotropy. 



and a minimum of 

e 



/NL 



-1.52- 



0.01 



2N ki0 (l + N kfi )+ 



+ ->/iV fc ,o(l + iV fc ,o)(l + 27V fc ,o) 



(23) 

where N = Nk,o- These values are depicted in Fig. Q]as 
a function of N k o- 



^ls = 14400 Mpc~ is the comoving distance to the sur- 
face of last scattering. We considered values of fc cut corre- 
sponding to a cutoff at I — 500 and / = 1500. We further 
calculated for 9 k ~ kr]o, as well as for 9 k ~ k/k* where 
k* <C k^ so that the terms with trigonometric functions 
in (flU)) average out as in (a) . The results for several val- 
ues of N k o are listed in Table H We see that the values 
for /nl are much smaller than for Nk ~ const. 



b) k cu t lies within the observable modes: 

As acknowledged earlier, if fc cut lies within the observ- 
able modes, there might be a conflict with measured 
power spectrum results; however, because of the impor- 
tance of fully investigating possible inflationary scenar- 
ios and understanding their consequences, we felt it was 
worthwhile to perform the calculation. Using a simple 
formula, one can relate the CMB-anisotropy multipole I 
to the wavenumber k most responsible for producing it 
at the surface of last scattering: k ~ l/ris [HI, where 



V. DISCUSSION AND CONCLUSION 

The conclusions from this work depend on the form of 
9k and the size of e. In the most conservative case, we 
suppose that 9k ~ r/ok and e < 0.01, and we find from 
(f2l"j) that /nl ^ 3 in the CMB and thus undetectable. 
This case conforms to the general idea that a detection of 
/nl would disfavor single-field inflation. However, if we 
allow either 9k to be constant or e to be slightly elevated 
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(which, as mentioned in Ei UVl is possible for a non-Bunch- 
Davies initial state), we can get a larger /nl that could 
be measured. The second of these scenarios - a larger e 
- can produce a (barely) measurable /nl ~ 15 if e ~ 0.1 
and Nk > 1. The first scenario, where 9 k ~ const, can 
produce a much larger and/or negative /nl, depending 
on the value of 9k, as we see from (f2"2")l . (j2"3")l . and Fig. 
ED rising to w 30 for 7V fc = 0(1) and e « 0.01. Thus, 
a detection of local-type /nl alone would not rule out 
single-find inflation. 

It might seem curious that / N l in ((2"Tj) docs not blow 
up for very large Nk t o- To understand this, note that 
(|4]) shows that /nl roughly measures the size of the 
bispectrum in the squeezed limit relative to the product 
of power spectra P{k\)P{k 3 ). However, for large Nk,o, 
P(k) oc Nk y o and Bq oc N% , so that Bq increases ex- 
actly like P{k\)P{k 3 ) with respect to Nk,o, and /nl is 
not enhanced. In fact, the non-numerical factors in (|21[) 
can be derived by supposing 



p(h)p(k 3 y 

In examining the effect of Nk.o of /nl in Table U we 
see that /nl decreases for decreasing / cut . This is consis- 
tent with expectations since the bispectrum we found was 
enhanced by k\/k 3 , and a lower l cut prevents the contri- 
bution of modes where k\ is large. Along the same lines, 
one might also wonder about the effect of other forms 
of Nk, for example a sharp peak at some momentum k 
due to a condensate or particle production. The effect of 
such a form would likely depend on the specifics of the 
model under consideration. If the boosted power is at 
small or large scales, where k\/k 3 is large, it will increase 
the observability of the /nl enhancement; however, the 
observability is also a function of the width of the peak, 
since an overly narrow peak will be washed out as the 



bispectrum is projected to 2D and then used to calculate 
/nl- 

As a last note, we remark on the fact that our results 
produce /nl larger than (5/12)(l — n s ) 0.01, as pre- 
dicted by the single-field consistency relation. To recon- 
cile this apparent discrepancy, we examine the integrals 
which give the enhanced terms, as described below (|13p . 
These integrals yield terms like e lfe3J, | jj . In the exact 
squcczcd-limit of the consistency relation, k 3 is taken to 
be arbitrarily small, so | A;3 770 1 1 and these terms are 
zero; the only pieces remaining reproduce the Maldacena 
squeezed- limit result. This explanation is somewhat un- 
satisfying because l^yol 1 implies that the initial con- 
dition was set when k 3 was far outside the horizon; how- 
ever, one may argue that unknown, early universe effects 
were responsible for this and, in any case, the consistency 
relation does hold. Thus, we see that the consistency re- 
lation is a powerful tool but it only applies perfectly in a 
certain limit (a similar conclusion was reached in [12j). 
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Appendix: Explicit Bispectrum 

From (|14|). we easily write out the bispectrum explic- 
itly: 



B/^(kx, fe, k 3 ) 
_i H 6 I 



(a kl + /3fcJ (a k2 + f3 k2 ) (a ks + /3 ks ) x 



v 



1 ff6 ^ / 1 \ 



4 02 fcifc 2 fc 3 

1 



ki + k 2 + k 3 
1 



k\ - k 2 + fc 3 
+ c.c. , 



k\ + k 2 - k 3 
1 



k\ - k 2 - k : 



(A.l) 
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where, in the last equation we let 77 — >• and averaged 
over r] , as discussed in § HT1 
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